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ABSTRACT 


Tie  NAVSTAR  Global  Positioning  System  is  a  worldwide, 
all-weather,  satellite  positioning  system  capable  of  high 
accuracy  real-time  position  determinations.  The  Applied 
Research  laboratories  (AfiL:UT),  in  conjunction  with  the 
Naval  Surface  Weapons  Center  (NSWC/DL)  and  other  government 
agencies,  conducted  geodetic  field  tests  of  a  government 
sponsored  prototype  receiver,  the  TI4100/GEOSTAR,  in  March 
1984.  Data  were  acquired  from  four  satellites  by  three 
receivers,  with  antennas  located  at  known  stations,  over 
approximately  six-hour  periods  each  day. 

Based  on  the  real-time  solutions  acquired  on  March  1,  8 
and  9,  1984,  absolute  point  position  determinations  have  an 
average  discrepancy  of  7.5  meters  with  a  one  sigma  repeat¬ 
ability  of  less  than  a  meter.  Relative  positions  were 
determined  to  an  average  relative  accuracy  of  1:9,700  for 
distances  of  14  to  26  kilometers. 
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I.  INTRODUCTION 


The  Naval  Surface  Weapons  Center  (NSWC) ,  with  sponsor¬ 
ship  ty  the  Defense  Happing  Agency  (DHA) ,  has  developed  a 
system  for  geodetic  positioning  with  tne  use  of  the  NAVSTAR 
Global  Positioning  System  (GPS)  £Ref.  1  ].  A  government 
sponsored  prototype  receiver,  the  TI4100/G2OSTAR,  was  devel¬ 
oped  by  Texas  Instruments  through  contract  with  Applied 
Research  laboratories.  The  University  of  Texas  at  Austin 
(A  RL :  OT)  [Ref.  2].  The  receiver  was  tasted  by  ARL:UT  in 
conjunction  with  NSWC  and  other  government  agencies  over 
short  and  medium  base  lines  in  March  1984.  The  Defense 
Mapping  Agency  Hydrographic/Topographic  Center  { DM AHTC)  has 
expressed  interest  in  independent  analyses  of  the  data  from 
these  tests  [Ref.  3]. 

Absolute  point  positioning,  as  it  is  used  in  this 
thesis,  is  the  determination  of  a  position  with  the  use  of 
one  GPS  receiver.  Relative  positioning  is  the  determination 
of  two  positions  relative  to  each  other  with  the  use  of  two 
receivers.  The  basic  premise  of  real-time  relative  posi¬ 
tioning  is  that  GPS  observation  errors  are  to  a  large  extent 
systematic  in  the  same  general  area  at  the  same  time.  If 
such  is  the  case,  relative  positioning  by  two  receivers  can 
locate  positions  relative  to  each  other  more  accurately  than 
absolute  positioning  alone. 

Accuracy  standards  for  horizontal  control  require 
certain  relative  accuracies  between  directly  connected  adja¬ 
cent  points.  Third-order  surveys  are  used  for  local  hori¬ 
zontal  control  in  hydrographic  surveying  and  other  local 
projects.  Third-order,  Class  II  standards  require  relative 
accuracies  of  at  least  1:5,000  between  directly  connected 
adjacent  points  (e.g.,  5  meters  relative  accuracy  between 
stations  separated  by  25  kilometers).  [Ref.  4] 


Geodetic  positioning  with  GPS  has  advantages  over 
conventional  surveying  methods.  Visibility  between  sites 
and  fair  weather  conditions  are  not  required.  Real-time 
positioning  has  a  special  application  in  the  case  where  the 
p^-jition  of  a  station  is  needed  in  a  short  amount  of  time 
and  before  the  data  can  be  post-processed.  An  analysis  of 
the  real-time  solutions  of  the  GEOSTAR  short  and  medium  base 
lice  tests  will  give  an  estimate  as  to  what  degree  of  accu¬ 
racy  the  system  can  achieve  in  real-time. 


II.  BACKGROUND 


A.  THE  GLOBAL  POSITIONING  SYSTEM 

The  Global  Positioning  System  is  a  worldwide,  all- 
weather,  satellite  positioning  system.  GPS  can  be  thought 
of  as  consisting  of  three  components;  the  Space  Segment,  the 
Control  Segment,  and  the  User  Segment.  The  Space  Segment 
will  consist  of  eighteen  satellites  (with  three  spares)  in 
six  orbit  planes  when  fully  operational  in  late  1S88  or 
early  1989  [Ref.  5].  The  satellites  will  orbit  the  Earth  in 
twelve-hour  periods  at  altitudes  of  approximately  20,000 
kilometers.  This  constellation  will  provide  four  satellites 
in  view  at  one  time  at  elevations  of  greater  than  20°  above 
the  horizon.  [Ref.  6]  The  satellites'  ground  paths  will  be 
fixed.  The  time  at  which  a  particular  satellite  will  be  in 
view  at  a  certain  position  will  be  about  four  minutes 
earlier  each  day  due  to  the  differences  between  solar  and 
sidereal  time.  [Ref.  7] 

The  GPS  satellites  have  on  orbit  weight  of  950  pounds. 
Power  is  supplied  by  solar  panels  and  nickel-cadmium 
batteries.  The  satellites  are  equipped  with  atomic  clocks 
with  stabilities  on  the  order  of  1  part  in  1.013.  [Ref.  8] 
Perturbations  in  the  sac01'ite  orbits  are  caused  by  the 
forces  of  attraction  of  astronomical  bodies,  gravitational 
irregularities  of  the  earth,  and  solar  radiation  pressure. 
These  perturbations  are  compensated  for  by  the  Control 
Segment  and  frequency  standard  offsets.  [Ref.  7] 

The  Control  Segment  consists  of  four  monitor  stations 
located  at  Alaska,  Hawaii,  Guam,  and  Vandenburg  AFB;  and  a 
Master  Control  Station  (MCS)  and  an  Upload  Station  at 
Vandenburg  AFB  [Ref.  9].  The  monitor  stations  track  the 


satellites  then  send  the  data  to  the  MCS  for  processing  to 
determine  updated  predictions  of  satellite  and  clock 
behavior.  The  updated  predictions  are  transmitted  to  the 
satellites  by  the  Upload  Station.  [Ref.  10] 

B.  TRANSMISSION  SIGNALS 

The  GPS  NAVSTAR  satellites  simultaneously  transmit  navi¬ 
gation  information  on  2  RF  frequencies;  the  LI  frequency  at 
1575.4  MHz.  and  the  L2  at  1227.6  MHz.  3oth  frequencies  are 
modulated  with  pseudo-random  noise  (PRN)  codes.  The  Precise 
Positioning  Service  (PPS)  ,  formerly  called  the  "P-code" ,  is 
modulated  on  both  frequencies  and  the  Standard  Positioning 
Service  (SPS) ,  formerly  called  the  "C/A-code",  is  available 
only  on  the  LI  frequency.  PRN  codes,  different  for  each 
satellite,  serve  the  purpose  of  uniquely  identifying  the 
particular  satellite  and  provide  a  means  of  determining  the 
signal  transit  time  by  measuring  the  phase  shift  necessary 
for  the  receiver  to  match  codes.  [ Eef .  11  and  12] 

The  SPS  is  a  relatively  short  code  of  1023  bits.  The 
code  repeats  itself  every  millisecond  at  a  bit  rate  of  1.023 
Mbps.  Due  to  the  short  duration  time  of  the  SPS,  the  code 
is  relatively  easy  to  match  and  lock  on  to.  The  PPS  has  a 
7-day  code  duration  and  10  times  the  bit  rate,  making  the 
PPS  code  relatively  difficult  to  match  with  the  receiver* s 
internally  generated  code.  The  PPS  offers  the  advantage  of 
a  finer  resolution  of  the  phase  shift  necessary  to  match 
codes,  thus  a  greater  measurement  acciv\cy  of  signal  transit 
time.  [Ref.  11] 

C.  THE  NAVIGATION  MESSAGE 

The  GPS  Navigation  Message  is  modulated  on  both  the  LI 
and  L2  frequencies  and  contains  information  necessary  to 
compute  the  user’s  position.  The  information  is  in  the  form 
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of  a  50  bps.  data  stream  common  to  both  tne  PPS  and  SPS. 

The  data  frame  is  1500  bits  long  and  consists  of  5  subframes 
of  6  seconds  (300  bits)  each.  The  Navigation  Message 
provides  information  such  as  satellite  e^hemerides,  system 
time,  satellite  clock  behavior,  and  transmitter  status. 

[Ref.  13] 

Subframes  are  divided  into  ten  30-bit  words.  The  first 
two  words  of  each  subframe  contain  system  time  and  time 
synchronization  information  for  transfer  from  the  SPS  to  the 
PPS  in  the  form  of  hand-over  words  (HOW)  and  telemetry  words 
(TLM) .  The  remaining  8  words  contain  user  navigation  data 
generated  by  the  Control  Segment.  [Ref.  13] 

Data  Block  1  occupies  words  2  through  10  in  subframe 
one.  It  contains  freguency  standard  corrections,  an  associ¬ 
ated  age  of  data  word,  and  ionospheric  delay  corrections  for 
single  freguency  receivers.  The  clock  correction  parameters 
are  used  to  determine  the  satellite's  clock  offset  from  GPS 
time  by  the  following  equations: 


t  =  t  -  At 

sv  sv 


(2.1) 


At 


sv 


=  a„ 


♦  an  (t  -  t  )  ♦  a.  (t  -  t  )  2 

1  oc  2  '  oc ' 


(2.2) 


where  t  is  the  GPS  time  in  seconds  measured  from  the  begin¬ 
ning  of  the  GPS  week  (approximately  Saturday/Sunday  midnight 
GMT),  t  is  the  satellite's  PRN  code  phase  time  at  the  tine 
of  message  transmission,  t  is  Data  Block  1  reference  time 
(seconds),  and  aQ,  a1#  and  a2  are  Data  Block  1  clock  correc¬ 
tion  parameters  (coefficients  of  a  second  order  polynomial)  . 
Note  that  t  may  be  approximated  by  t  in  equation  2.2. 

[Ref.  13] 

The  age  of  data  word  (AODC)  is  the  time  between  Data 
Block  1  reference  time  and  the  time  since  the  last  measure¬ 
ment  used  in  computing  the  clock  correction  parameters 


[Ref.  13].  Test  data  in  this  thesis  ware  obtained  using 
both  the  LI  and  L2  frequencies,  therefore,  the  ionospheric 
delay  parameters  for  single  frequency  receivers  contained  in 
Data  Block  1  are  not  used  (see  Van  Dierendonck,  et  al, 

1978)  . 

Data  Block  2  contains  ephemeris  information  in  subframes 
2  and  3.  This  data  includes  associated  age  of  data  words 
(A3DE)  which  represent  the  time  difference  between  Data 
Block  2  reference  time  and  the  time  of  the  last  measurement 
used  in  the  estimation  of  the  ephemeris  parameters.  The 
ephemeris  data  includes  conventional  Keplerian-t ype  parame¬ 
ters,  secular  drift  terms,  and  harmonic  coefficients  used  to 
cocrect  for  perturbations  in  the  Kaplerian  orbit.  [Ref.  13] 
The  satellite's  position  is  calculated  by  Van  Dierendonck ' s 
algorithm  [fief.  13]  with  the  addition  of  a  secular  drift 
term  to  correct  for  perturbations  in  the  inclination  of  the 
satellite's  orbit  [Ref.  14]. 

The  Message  Block  appears  in  the  third  through  the  tenth 
30-bit  words  of  subframe  4  of  the  navigation  Message.  This 
space  is  reserved  for  alphanumeric  messages  provided  to  the 
user  by  the  Control  Segment.  Space  exists  for  23  eight-bit 
ASCII  characters.  [Ref.  13] 

Almanac  information  appears  in  Data  Block  3  in  the  fifth 
subframe.  Data  Block  3  consists  of  25  subframes  allowing 
for  almanacs  of  up  to  24  satellites  (with  one  dummy 
subframe).  Twenty-five  navigation  messages  must  be  received 
to  acquire  the  complete  data  set  of  Data  Block  3.  Extra 
subframes  can  be  used  to  repeat  the  almanacs  for  certain 
satellites.  Satellite  almanacs  provide  the  user  with  less 
precise  position  and  clock  correction  information.  Almanacs 
aid  the  user  in  satellite  selection  and  simplify  the  acqui¬ 
sition  of  signals.  Almanacs  include  ephemeris  parameters, 
clock  correction  parameters,  satellite  identification,  and 
satellite  health  terms.  [Ref.  13] 


D.  RECEIVES  OPERATIC! 


The  GEOSTAR  GPS  receiver  is  a  government  sponsored 
version  of  the  TI4100  commercial  receiver  (Fig.  2.1).  The 
GEOSTAR  is  a  single  channel,  multiplex  receiver  capable  of 
tracking  four  satellites  on  the  11  and  L2  PPS.  The  control 
and  display  unit  (CDO)  is  a  hand-held  input/output  device 
with  two  16-character  alphanumeric  displays.  The  antenna 
assembly  consists  of  a  conical  omnidirectional  antenna  with 
a  preamp  that  allows  remote  operation  of  the  receiver.  Two 
100-foot  cables  connect  the  antenna  assembly  to  the 
receiver.  One  cable  is  used  to  input  DC  power  to  the  preamp 
and  output  the  amplified  RF  signals.  The  other  cable  is 
used  for  a  built-in  test  signal  injection  from  the  receiver. 
[Ref.  15] 

A  typical  GEOSTAR  geodetic  surveying  unit  consists  of  a 
2-man  team,  a  station  wagon  or  van,  and  the  GEOSTAS  field 
set.  The  field  set  contains  the  receiver,  CDU,  the  antenna 
assembly,  a  cassette  tape  recorder,  a  tripod  with  a  leveling 
bracket,  two  12-volt  automobile  batteries,  an  EFRATOM  FRK-H 
rubidium  oscillator,  a  weather  station  with  sensors,  cables, 
transportation  cases,  tools,  manuals,  and  spare  parts.  The 
field  set  is  easily  contained  and  transported  in  the  rear  of 
the  vehicle. 

The  antenna  assembly  is  mounted  on  a  standard  surveyor's 
tripod  with  a  leveling  bracket  then  leveled  and  plumbed  over 
the  mark  to  be  surveyed.  A  reference  mark  on  the  base  of 
the  antenna  assembly  is  oriented  to  north  using  a  hand 
compass.  The  height  of  the  L2  phase  center  mark  at  the  base 
of  the  antenna  assembly  above  the  survey  point  is  measured 
and  recorded.  The  antenna  assembly,  the  power  supply,  tne 
cassette  recorder,  the  weather  station,  and  the  rubidium 
oscillator  are  then  connected  to  the  receiver  with  the 
appropriate  cables.  Scratch  tapes  are  inserted  in  the 
cassette  recorder  for  the  power-up  procedure. 
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The  system  is  turned  on  and  the  receiver  and  cassette 
recorder  self-tests  are  initiated.  When  the  self-tests  are 
successfully  completed,  the  GESAR  bootstrap  tape  is  inserted 
in  Drive  1.  The  program  is  loaded  and  blank  tapes  are 
inserted  in  the  two  tape  drives.  The  operator  is  then 
prompted  by  the  CDU  for  input  information  including  an  esti¬ 
mate  of  the  user's  position,  weather  data,  satellites  to  be 
tracked,  and  Doppler  aiding  values  for  each  satellite.  The 
system  then  operates  unattended  except  for  the  replacement 
of  cassette  tapes  as  required  and  the  periodic  manual  input 
of  weather  data  (if  the  automatic  weather  station  is  not 
used)  . 

E.  POSITION  CONFUTATION 

The  satellite's  position  and  the  time  of  signal  trans¬ 
mission  are  computed  with  the  broadcast  ephemeris  and  the 
clock  correction  parameters  included  in  the  Navigation 
Kessage.  The  time  between  the  transmission  and  the  recep¬ 
tion  of  the  satellite  signal  multiplied  by  the  speed  of 
lignt  is  the  pseudorange. 

Pseudoranges  used  to  calculate  the  real-time  solutions 
are  corrected  for  atmospheric  delays  in  the  GEOSIAR  Receiver 
[fief.  17].  Ionospheric  refraction  is  nearly  inversely 
proportional  to  the  square  of  the  frequency.  By  receiving 
both  the  Li  and  the  L2  frequencies,  the  ionospheric  delay  is 
calculated  by  comparison.  [Ref.  16]  Tropospheric  delay  is 
calculated  by  a  model  based  on  temperature,  pressure, 
humidity,  and  the  elevation  of  the  satellite  above  the 
horizon  [fief.  17].  Knowing  the  positions  of  four  satellites 
and  the  ranges  to  the  user,  the  user's  position  is  calcu¬ 
lated  by  a  least  squares  method. 


F.  GEODETIC  DATUMS 


The  GEOSTAR  real-time  solutions  are  referenced  in  DoD  World 
Geodetic  System  1972  (WGS-72)  Cartesian  coordinates.  WGS-72 
is  an  earth-centered,  earth-fixed,  glonal,  geodetic  datum 
adopted  by  DMA  in  March  1974  [Ref.  18].  The  reference 
surface  is  an  ellipsoid  that  best  approximates  the  shape  of 
the  earth’s  geoid.  The  ellipsoid  is  normally  defined  by  the 
length  of  the  semimajor  axis  and  the  flattening. 

Latitude  ranges  between  -90°  and  90°  with  north  posi¬ 
tive.  Longitude  is  measured  positive  east  from  0°  to  just 
less  than  360°  from  the  WGS-72  reference  meridian  at  0".54 
east  of  the  Greenwich  Meridian.  The  Cartesian  X-axis  is 
formed  by  the  intersection  of  the  reference  meridian  plane 
and  the  plane  of  the  equator.  The  Y-axis  is  90°  east  of  the 
X-axis  in  the  plane  of  the  equator.  The  Z-axis  is  in  the 
direction  of  the  earth's  mean  rotation  axis  based  on  the 
average  terrestrial  pole  of  1900-05.  [Ref.  19] 

The  North  American  Datum  1927  (NAD-27)  is  a  regional 
datum  based  on  the  Clarke  1866  Ellipsoid  with  the  origin  at 
Meades'  Ranch,  Kansas.  The  Clarke  1866  Ellipsoid  minimizes 
tne  deflection  of  the  vertical  and  geoidal  undulations  in 
the  continental  United  States.  The  National  Ocean  Survey 
(NOS) ,  the  National  Geodetic  Survey  (NGS) ,  and  the  United 
States  Geological  Survey  (U5GS)  use  NAD-27  coordinates  in 
most  survey  operations.  The  transformation  from  NAD-27  to 
KGS- 7 2  reference  systems  creates  offsets  of  the  order  of  5 
meters  [Ref.  2].  It  is  assumed  that  this  error  is  removed 
when  the  reverse  transformation  is  made;  but  the  error 
source  will  still  exist  in  the  analysis  of  the  WGS-72  point 
position  discrepancies. 


III.  AHALYSIS  METHODS 


A.  EIPEBIMENT  DESIGH 

Relative  position  field  testing  of  tne  3SOSTAR  Receiver 
was  performed  in  central  Texas  in  February  and  larch  193a. 
Three  receivers,  each  with  an  antenna  located  over  a 
different  known  position,  were  operated  continuously  over 
approximately  the  same  period  of  satellite  visibility.  Four 
satellites,  with  PEN  codes  6,  3,  9,  and  11,  were  tracked  in 
these  tests.  The  period  of  data  gathering  was  coordinated 
with  the  time  period  in  which  all  four  satellites  were  in 
view  at  the  same  time.  The  real-time  positions  of  the 
antennas  were  recorded  every  12  seconds  over  approximately  a 
six-hour  period  each  day.  The  data  were  recorded  on 
cassette  tape  then  transferred  to  9-track  tape  by  ARL :  (JT. 

B.  GPS  USER  TAPES 

The  recorded  data  and  the  associated  position  coordi¬ 
nates  used  in  the  analysis  were  obtained  from  NSUC/DL.  Tne 
GPS  User  Tapes  are  a  standardized  format,  9-track  magnetic 
tape.  The  tapes  are  unlabelled,  ASCII  coded,  1600  bpi,  with 
80  characters  per  physical  record  and  1  record  per  block. 

The  first  file  on  the  tapes  consists  of  header  comment 
records.  This  file  assists  the  user  in  identifying  the  tape 
and  includes  information  such  as  the  date  and  site  of  the 
data  acquisition. 

Data  files  begin  with  an  80-character  control  record 
(Data  Item  List)  which  identifies  the  types  of  data  tnat 
follow.  The  first  location  of  all  records  is  reserved  for  a 
control  character.  The  Data  Item  List  controL  record  has  an 
asterisk  **'  in  the  first  position  followed  by  an  integer  in 


the  next  three  positions  which  specifies  the  numcer  of  times 
to  repeat  the  sequence  of  the  Data  Item  List.  The  ren.airu.n-j 
76  characters  contain  up  to  IS,  4-digit  data  type  identi¬ 
fiers  which  specify  the  type  of  data  that  is  to  follow  the 
control  record.  Control  records  are  read  with  a 
(A  1 , 13,1914)  format. 

Data  information  are  one  or  more  records  in  length  and 
divided  into  groupings  called  data  items  (Tab.  I) 

[Ref.  19]. 

TABLE  I 

Data  Types  of  GPS  User  Tapes 
Item  Identifier  Data  Type 

0001-0999  Exchange  conversion  program  information 

1000-1999  Constant  equipment  related  data 

2000-2999  Campaign  identifying  data 

3000-3999  Infrequent  measurement  data 

4000-4999  Space  vehicle  related  data 

5000-5999  Weather  data 

6000-6999  Real-time  solution  lata 

7000-7999  Measurement  data 

8000-8999  Post  analysis  data 

The  data  type  identifiers  of  the  Data  Item  List  specify  the 
order  and  the  format  of  how  the  data  is  to  be  read. 

Data  used  in  the  analysis  were  read  from  the  user 
exchange  tapes  using  a  slight  modification  of  the  sample 
input  program  contained  in  reference  19  (App.  A). 

Real-time  solution  information  (data  type  6010)  was  read 
from  the  tapes  and  written  to  mass  storage  and  retained. 
Other  information  such  as  measurement  data  types  (data  type 
920),  real-time  solution  configuration  (data  type  1150),  and 


antenna  location  (data  type  2210)  were  read  manually  from 
the  first  page  of  the  printouts  of  the  tapes. 

C.  COORDINATE  TRANSFORMATIONS 

The  coordinates  of  the  stations  used  in  the  field 
testing  were  given  in  NAD-27  geodetic  coordinates.  The 
real-time  solution  output  of  the  GEOSTAfi  Receiver  are  in 
WGS-72  Cartesian  coordinates.  To  make  comparisons  of  the 
given  positions  to  the  observed  positions,  the  given  posi¬ 
tions  were  transformed  to  WGS-72  Cartesian  coordinates.  The 
transformation  was  done  using  the  Abridged  Mol  indensky 
formulas  contained  in  the  NGS's  conversion  program  (D035-41) 
for  the  HP-4 1CV  programmable  calculator.  The  program 
prompts  the  user  for  the  values  of  the  origin  shift  between 
coordinate  systems.  The  X,  Y,  and  Z  origin  shifts  for 
transforming  from  NAD-27  to  WGS-72  are  -22,  157,  and  176 
meters,  respectively  [Ref.  20].  These  values  differ  from 
the  origin  shift  values  supplied  with  the  transformation 
program  by  NGS.  The  NGS  values  are  -12.547,  155.804,  and 
175.369,  for  the  X,  Y,  and  Z  origin  shifts,  respectively. 

The  program  generates  semimajor  axis  and  inverse  flattening 
parameters  of  both  reference  systems  (Tat.  II) . 


TABLE  II 

Transf oraation  Program  Parameters 
Ellipsoid  a  (meters)  1/f 


Clarke  1866 
WGS  72 


6,378,206.4 

6,378,135.0 


294. 9786982 
298.2600000 
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D.  SATELLITE  UPDATES 

The  ephemeris  elements  and  clock.  model  updates  occurred 
at  identical  times  according  to  the  data  recorded  in 
field  tests.  An  approximation  of  the  update  times  is  calcu¬ 
lated  by  subtracting  the  Age  of  Data  (AODE)  from  tae  ephem- 
eris  reference  time  (t  ) ,  i.e., 

update  time  =  t  -  AODE  (3.1) 

oe 

The  update  times  changed  by  -4  96  ,  1052,  or  1552  seconds  when 
the  computation  of  equation  3. 1  was  performed  for  each 
satellite  with  subsequent  occurrences  of  the  broadcast 
ephemeris  on  the  recording.  The  change  is  due  to  packing  of 
bits  in  the  transmissions  of  the  AODE  values.  An  approxima¬ 
tion  of  the  update  times  is  made,  within  a  few  minutes,  by 
using  the  first  AODE  values  that  cnange  to  smaller  numbers 
(on  the  order  of  1  to  2  hours)  from  previously  recorded 
large  values  (on  the  order  of  24  hours).  [Bef.  14] 

E.  A VEB AGING  HETHODS 

Relative  position  determinations  were  made  using  real¬ 
time  solutions  from  two  different  stations  at  the  same  GPS 
time  tags.  Computer  programs  were  written  to  read  the  real¬ 
time  solutions  from  mass  storage  and  compute  the  discrepan¬ 
cies  between  the  observed  and  the  given  WGS-72  coordinates 
of  pairs  of  stations  at  common  GPS  times.  The  beginning  and 
ending  times  of  data  acquisition  varied  among  the  three 
receivers  operated  during  the  same  time  period;  thus,  the 
effect  of  comparing  solutions  at  the  same  time  tags  resulted 
in  the  elimination  of  some  of  the  real-time  solutions  from 
the  analysis.  In  effect,  the  tape  with  the  earlier  starting 
time  was  spooled  forward  to  correspond  to  the  starting  time 
of  the  first  real-time  solution  recorded  by  the  receiver 
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that  began  operation  at  a  slightly  later  time.  Similarly, 
the  time  of  the  last  solution  common  to  both  receivers  as 
tne  last  solution  recorded  by  the  receiver  that  ended  opera¬ 
tion  first. 

A  mean  position  was  determined  in  two  different  ways  to 
show  the  effect  of  the  ephemeris  and  clock  updates  on  the 
real-time  solutions  and  to  eliminate  real-time  solutions 
with  large  discrepancies  from  the  averaging.  In  general, 
relatively  large  discrepancies  are  evident  in  roughly  the 
first  hour  of  receiver  operation.  The  real-time  solution 
discrepancies  of  pairs  of  stations  were  averaged  beginning 
with  the  values  one  hour  after  and  two  hours  after  the  first 
solution  time  common  to  both  receivers.  In  these  tests, 
beginning  the  averaging  one  hour  after  the  initial  solutions 
corresponds  to  averaging  the  solutions  that  were  recorded  a 
short  time  before  or  after  the  update  times  and  through  the 
remainder  of  the  recording  session.  Beginning  the  averaging 
two  hours  after  the  initial  solutions  corresponds  to  using 
only  those  values  that  were  recorded  much  after  the  updates 
and  through  the  remainder  of  the  session.  Real-time  solu¬ 
tions  are  treated  as  independent  observations  in  the  statis¬ 
tical  analysis. 

F.  COMPUTATIONS 

The  distances  between  stations  (d)  ,  in  the  W3S-72  refer¬ 
ence  system,  are  computed  by  eguation  3.2  where  XL  ,Ylf  and 
Z1  are  the  RGS-72  Cartesian  coordinates  of  one  of  the 
stations,  and  X9 ,  Y 2 ,  and  Z 9  are  the  coordinates  of  the 
other  station  of  the  side.  The  side-lengths  are  calculated 
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for  both  the  given  coordinates  and  the  mean  values  of  the 
GEOSTAE  real-time  solution  coordinates.  The  discrepancies 
between  the  two  are  computed  by  subtracting  the  given 
distances  from  the  observed  distances.  Discrepancies  are 
used  in  the  averaging  to  save  computer  time  and  to  reduce 
round-off  error. 

The  absolute  value  of  the  point  position  discr epancies 
are  computed  by  equation  3.3  where  the  subscripts  o  and  g 
represent  the  observed  and  given  values,  respectively. 

5  =  y<x0  -  v2  *  <yo  -  v2  +  (zo  -  v2  <3-3) 

The  standard  error  of  a  single  observation  was  calcu¬ 
lated  as  a  measure  of  the  precision  of  the  real-time  solu¬ 
tions.  The  means  and  standard  errors  of  a  single,  12-second 
real-time  solution  were  calculated  for  the  coordinate 
discrepancies  and  the  point  position  discrepancies. 
Comparisons  are  made  for  the  two  averaging  methods  (Sect. 

Ill  E) .  The  mean  WGS-72  Cartesian  coordinate  discrepancies 
were  used  to  compute  the  final  observed  positions.  The 
observed  positions  are  computed  by  subtracting  the  discrep¬ 
ancies  from  the  given  coordinates. 

G.  GEODETIC  COMPARISONS 

Geodetic  accuracy  standards  for  horizontal  control 
require  certain  relative  accuracies  between  directly 
connected  adjacent  points.  The  distances  between  points  are 
computed  along  the  ellipsoid  -  the  geodesic.  To  compare  the 
accuracies  of  the  observed  positions  with  the  standards  for 
geodetic  control,  the  WGS-72  Cartesian  coordinates  of  the 
observed  positions  are  transformed  to  NAD-27  geodetic  coor¬ 
dinates.  The  transformations  were  done  using  a  modification 
of  NGS’s  D035  program  for  the  HP-4 1 C V  programmable 
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calculator.  The  modified  version  uses  the  same  parameters 
as  those  contained  in  the  original  program  (Tab.  II). 

The  geodetic  inverse  computation  is  used  to  calculate 
the  geodetic  distance  and  azimuths  between  two  known  points. 
The  inverses  between  given  positions  and  between  observed 
positions  were  computed  and  compared. 

To  simulate  the  positioning  of  an  unknown  station,  one 
station  of  a  side  is  held  fixed  and  the  other  station  is 
located  using  the  direct  computation.  The  coordinates  of 
the  known  position  and  the  distance  and  azimuth  to  the 
unknown  position  are  the  input  variables  and  the  coordinates 
of  the  unknown  station  are  computed.  The  distance  and 
azimuth  calculated  by  the  inverse  computation  of  the 
observed  point  positions  are  used  in  tnis  calculation.  The 
direct  and  inverse  computations  were  done  on  the  IBM  3033 
based  on  Puissant's  Coast  and  Geodetic  Survey  Formula 
[Ref.  21  ]. 

The  standard  errors  of  the  discrepancies  of  the  geodetic 
position  determinations  for  entire  recording  sessions  are 
calculated  as  an  approximation  of  the  repeatability.  To  do 
so,  the  two  point  position  discrepancies  of  each  station  on 
days  61  and  63  must  first  be  averaged  to  determine  a  single 
value  for  the  station. 
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IV.  BESPITS  AND  DISCPSSION 


A.  GIVEN  POSITIONS  AND  DISTANCES 

Data  in  this  report  are  from  relative  position  tests  of 
March  1#  8,  and  9,  1984.  Three  receivers  were  operated 
simultaneously  over  a  period  of  approximately  six  hours. 

The  antennas  of  the  March  1st  tests  were  located  over  a 
second-order  position  and  two  of  its  reference  marks.  The 
triangle  has  side- lengths  of  approximately  30  meters.  The 
tests  of  March  8th  and  9th  were  conducted  at  first-order 
survey  points  separated  by  distances  of  14  to  25  kilometers. 
The  data  of  the  tests  at  one  of  the  positions  on  March  9th 
was  not  available  at  the  time  of  this  writing.  The  data 
gathered  on  March  1st  will  be  defined  as  the  short  base  line 
tests  and  the  14  to  25  kilometer  testing  will  be  termed 
medium  base  line  tests.  Tables  III  and  IV  give  the  names 
and  the  NAD-27  coordinates  of  the  stations  and  Table  V  gives 
the  distances  and  azimuths  in  each  triangle.  Distances  and 
azimuths  were  calculated  using  the  inverse  computation 
(Sect.  Ill  G)  . 

The  NAD-27  coordinates  of  the  given  positions  were 
transformed  to  KGS- 72  Cartesian  coordinates  using  a  NGS 
transformation  program  (D035-4  1).  The  ellipsoidal  heights 
are  corrected  for  the  heights  of  the  antennas  above  the 
station  marks  prior  to  the  transformations.  The  heights  of 
the  antennas  above  the  stations  (data  type  2110)  were  read 
manually  from  the  first  page  of  the  printouts  of  the  tapes. 
The  distances  between  the  points  in  W3S-72  Cartesian  coordi¬ 
nates  were  calculated  using  equation  3.2  (Tab.  VI). 

The  computed  side  distances  for  the  FGS-72  Cartesian 
coordinate  system  are  in  all  cases  greater  than  or  equal  to 


the  computed  geodetic  distances.  In  the  medium  base  line 
tests,  the  differences  were  found  to  increase  by  .0035 


TABLE  III 

Station  Names 

J.D. 

Number 

Name 

Abbre via  tion 

61 

85019 

ARL  10563  DMA 

JME 

35020 

ARL  10563  RM  1  DMA 

RM  1 

35021 

ARL  10563  RM  2  DMA 

RM2 

68,69 

35022 

BF.  ZERO 

3ZERG 

85027 

PILOT  KNOE  1935 

PILOT 

*  85028  BALD  KNOB  1935  BNOB 


*  Station  not  used  in  March  9th  analysis 


TABLE  IT 

NAD-27  Station  Coordinates 


St  ation 

North  Latitude 

West 

Longitude 

H  t .  ( m . ) 

JMR 

30° 

22’ 

56". 74285 

970 

43’ 

35". 97558 

235.  5 

EMI 

30  0 

22’ 

57". 54452 

970 

43’ 

36". 31041 

235.  7 

FM  2 

30° 

22' 

57". 06234 

970 

43’ 

35". 01009 

235.6 

BZERC 

30° 

23’ 

16". 31238 

970 

43’ 

42". 05790 

238.  4 

FILOT 

300 

OS' 

30". 93387 

970 

42’ 

28". 08922 

215.9 

END  3 

300 

20' 

07". 61 000 

970 

35’ 

30". 84700 

204.  4 

percent 

of  the 

length  of  the 

line. 

Differences 

of  .911, 

.748,  and  .534  meters  were  computed  for  the  side-lengths  of 
approximately  25.5,  22.6,  and  14.3  kilometers,  respectively. 
These  differences  occur  because  the  Cartesian  distances  are 
computed  between  each  point;  the  geodetic  distances  are 
computed  from  the  projections  of  the  points  on  the  ellipsoid 
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which  does  r.ot  accour.t  for  the  ellipsoidal  heights.  The 
distances  between  BZERO  and  PILOT  are  equal  to  the  nearest 


TABLE  7 

NAD-27  Geodetic  Distances  and  Azimuths 


Side 

Distance {m. ) 

Azimuth 

JMR  -  RSI 

26.255 

340° 

05’ 

34". 72573 

J ME  -  RM2 

27. 591 

69° 

06’ 

37"  .  54898 

EMI  -  RK2 

37.759 

113° 

09' 

1 9"  .  14446 

BZERO  -  PILOT 

25,492.277 

1750 

32’ 

48". 25375 

PZERO  -  BNOB 

14,346. 955 

1130 

51* 

26" . 80840 

BNGP  -  PILOT 

22,556. 303 

209° 

40 » 

09". 53915 

TABLE  VI 


WGS-72  Side-length 

Distances 

J.  0 .  Side 

Distance (a. ) 

61  JMR  -  P, Ml 

26.255 

JMR  -  R M2 

27.593 

RM1  -  RM2 

37.759 

63  BZEPO  -  PILOT 

25, 433 . 188 

BZERO  -  B NOB 

14, 347.439 

BNOB  -  PILOT 

22, 557.051 

69  BZERO  -  PILOT 

25,493.188 

millimeter 

for  tests  run  on  days  69 

and  59  even 

though  the 

heights  of 

the  antennas  above  the  s 

tations  were 

different  on 

both  days. 
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B.  SATELLITE  UPDATE  TIMES 

The  ephemeris  and  clock  update  tines  occurred  approxi¬ 
mately  one  hour  after  the  start  of  data  acquisition  in  tnese 
tests.  There  are  approximately  25  to  30  recordings  of  the 
broadcast  ephemeris  (data  type  4110)  during  the  six-hour 
recording  periods.  The  update  times  were  computed  (eqn. 

3.1)  using  the  AODE  values  that  changed  from  large  numbers 
to  much  smaller  numbers  in  subsequent  recordings  of  the 
broadcast  ephemeris  (Tab.  VII).  Satellites  with  PRN  codes 
9  and  11  were  updated  approximately  one-half  hour  after  the 
updates  of  codes  6  and  8  on  days  61  and  69.  According  to 
the  recorded  data,  all  four  satellites  were  updated  at 
approximately  the  same  time  on  day  68. 

TABLE  VII 

Satellite  Update  Times  (Seconds) 


J.  D. 

PFN  codes  t 

AODE 

t_.  -  AO 

oe  ut 

61 

6,8 

378,000 

6,  144 

371, 856 

9,11 

378,000 

4,  096 

373,904 

68 

6,9,9,11 

378,000 

6,144 

371,856 

69 

e  P 

'  t 

460,300 

4,  096 

456, 704 

9,11 

464,400 

6,  144 

458,256 

C.  BEHAVIOF.  DF  THE  REAL-TIME  SOLUTIONS 

1 .  WGS-72  Coordinate  Discrepancies 

The  real-time  solutions  were  read  from  the  GPS  User 
Tapes  and  retained.  The  solutions  of  pairs  of  stations, 
with  time  tags  of  real-time  solutions  common  to  both,  were 


» 
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used  in  the  analysis.  The  discrepancies  between  the 
observed  coordinates  and  the  given  W3S-72  coordinates 
{observed  -  given)  were  computed  and  plotted  over  time 
(Figs.  4.1  -  4.8).  The  vertical  arrows  in  the  figures 
indicate  the  times  of  satellite  updates  (Tab.  VII).  The 
tick  marks  on  the  horizontal  (time)  axis  represent  one  hour. 

In  general,  the  X,  Y,  and  Z  coordinate  discrepancies 
begin  with  large  values  in  the  initial  solutions  and  then 
become  smaller  and  more  consistent  after  satellite  updates 
and  over  time.  There  is  a  significant  improvement  in  the 
discrepancies  of  the  observed  coordinates  near  the  time  of 
the  update  of  PRN  codes  6  and  8  at  station  JdR  and  near  the 
update  of  codes  9  and  11  at  station  P'll  on  day  61  (Figs. 

4.1  and  4.2).  The  precision  of  the  discrepancies  in  each 
coordinate  increase  at  ail  three  stations  on  day  61,  a  snort 
time  after  the  update  of  codes  9  and  11.  The  discrepancies 
at  stations  BZEEO  and  PILOT  on  day  68  decrease  approximately 
one  hour  after  the  initial  solutions  (Figs.  4.4  and  4.5). 

2.  WGS-72  Point  Position  Discrepancies 

The  absolute  values  of  the  point  position  discrepan¬ 
cies  of  pairs  of  stations  with  common  time  tags  were 
computed  (ecn.  3.3)  for  five-minute  time  segments  and 
plotted  with  the  standard  errors  over  time  (Figs.  4.9  - 
4.16).  In  general,  the  point  position  discrepancies  and  the 
standard  errors  decrease  after  satellite  updates  and  over 
time.  Except  for  the  general  trend,  there  is  no  obvious 
correlation  of  the  accuracies  (discrepancies)  and  the  preci¬ 
sion  of  the  five-minute  time  segments  at  all  stations  on  the 
same  day. 

3*  WGS-72  Side-length  Discrepancies 

Felatively  high  accuracy  side-length  determinations 
are  evident  at  certain  times;  however,  the  times  of  these 


30 


DISCREPANCY  (METERS) 

-70.0-60  0  -50.0  -  40.0  -  30.0  -20.0  -10.0  0.0  10.0  20.0  30.0  4  0.0  50.0  60.0  70.0 


ME  I  ERS) 

10.0  20.0  30.0  40.0  50.0  60.0  70.0 


L(j.  ' '90-1 

MEDIUM  BADE  LINE  "ErT 
SAMPLING  INTERVAL  =  5  MING 'EE 


Figure  4.4  HGS-72  Coordinate  Discrepancie 
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Figure  4.6 


fGS-72  Coordinate  Discrepancies 
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J.D.  61,1984 

SHOPT  BASELINE  TES^S 
5  MINUTE  SEGMENTS 


Figure  4.12  Discrepancies  and  Standard  Errors 
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STATION  PILOT 
J.D.  68,I98T 


Figure  4.13  Discrepancies  and  Standard  E 
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Figure  4.14  Discrepancies  and  Standard  Errors 
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solutions  are  not  consistent  at  all  stations  luring  the  same 
day  (Figs.  4.17  -  4.23).  The  discrepancies  in  the  lengths 
of  sides  JMF.-PM2  and  RM1-RM2  are  relatively  snail  near  time 
tag  282632;  however,  the  side- length  iiscrepancy  of  JMR-RS1 
at  approximately  the  same  time  is  relatively  larger  than 
other  determinations  of  that  side  on  lay  61.  The  side- 
length  discrepancies  of  JMP-RM2  and  RM1-EM2  are  -.044  ana 
.037  meters  at  time  tags  332632  and  332600,  respectively. 

The  discrepancy  of  JfiP-PMl  at  time  tag  3S2600  is  .955 
meters;  whereas,  the  discrepancy  is  .476  meters  at  time  tag 
386400.  Therefore,  the  most  accurate  side-length  determina¬ 
tions  do  not  necessarily  occur  at  approximately  the  same 
time  for  all  three  sides. 

The  most  accurate  side-length  determinations  do  not 
necessarily  occur  when  the  magnitudes  of  the  point  position 
discrepancy  vectors  of  two  stations  are  nearly  the  same 
because  of  differences  in  the  directions  of  the  discrepan¬ 
cies.  The  side-length  discrepancies  determined  on  day  66 
for  the  side  8ZEF.C-8NOB  illustrate  this  point  (Fig.  4.21). 
The  point  position  discrepancies  for  stations  3ZERO  and  3N0B 
are  10.889  and  10.816  meters  at  time  tag  377808  and  the 
side-length  discrepancy  is  -.151  meters.  At  time  tag 
384408,  the  side-lenyth  discrepancy  is  -.036  meters  and  the 
magnitudes  of  the  point  position  discrepancies  are  9.295  and 
14.053  meters.  The  side-length  discrepancy  is  smaller  at 
the  later  time  tag  even  though  the  difference  in  the  magni¬ 
tudes  of  the  point  position  discrepancies  are  greater  at 
that  time. 

D.  AVERAGING  RESULTS 

1 •  Determination  of  the  Method  of  Averaging 

Two  averaging  schemes  were  used  to  arrive  at  mean 
position  discrepancies  (Sect.  Ill  E) .  The  means  and 
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OLUTE  VALUE  OF  DISCREPANCY  (.METERS) 


Figure  4.23  Point  Position  and  Side-length  Discrepancies 
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standard  errors  of  discrepancies  were  computed  for  the  X,  Y, 
and  1  coordinates  and  the  absolute  values  of  the  total  point 
position  discrepancies  (Tabs.  VIII  and  IX). 

The  standard  errors  of  the  point  position  discrepan¬ 
cies  for  the  solutions  that  were  averaged  beginning  two 
hours  after  the  initial  solutions  are,  in  general,  smaller 
than  those  that  were  determined  using  all  solutions  begin¬ 
ning  one  hour  after  the  initial  solutions.  The  average 
standard  error  of  the  absolute  value  of  the  point  position 
discrepancies  for  all  pairs  of  determinations  is  2. 302 
meters  using  the  former  averaging  method,  compared  to  1.422 
meters  using  the  latter  averaging  method.  Much  of  the 
difference  between  these  two  standard  errors  can  be  attrib¬ 
uted  to  the  relatively  large  point  position  discrepancies 
computed  for  stations  3ZEEO  and  PILOT  on  day  69  between  the 
ficst  and  second  hours  after  the  initial  solutions  (Figs. 
4.15  and  4.16).  Neglecting  stations  BZERO  and  PILOT  on  day 
60,  the  average  of  all  standard  errors  of  the  point  position 
aiscr epancies  are  1.657  meters  and  1.333  meters  for  the  two 
different  averaging  methods. 

The  elimination  of  the  real-time  solutions  recorded 
between  one  hour  and  two  hours  after  the  start  of  data 
acquisition  results  in  smaller  standard  errors  of  the  point 
position  discrepancies  in  the  averaging  process.  The 
remainder  of  the  analysis  uses  the  mean  coordinates  deter¬ 
mined  by  the  averaging  method  in  which  only  the  real-time 
solutions  recorded  two  hours  after  the  initial  solutions  are 
used  (Tab.  IX)  . 

2 .  WGS-72  Results 

The  Y  coordinate  discrepancies  are  larger  than  the 
discrepancies  of  X  and  Z  coordinates  at  all  stations  with 
the  exception  of  those  at  station  BZERD  on  day  59.  The  mean 
observed  X  coordinates  are  consistently  different  from  the 
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TABLE  IX 

Means  and  Sigmas  Beginning  Two  Hours  After  Initial 
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giver.  X  coordinates  Ly  approximately  5  to  6  meters  at  all 
stations  on  all  days  of  data  acquisition  Laving  an  average 
standard  error  of  .564  meters  for  all  determinations.  The  Z 
coordinate  has  smaller  discrepancies  at  ail  stations  on  days 
68  and  69  compared  to  the  other  coordinates.  The  average 
standard  errors  of  the  Y  and  Z  coordinate  discrepancies  for 
all  determinations  are  1.600  and  1.560  meters,  respectively 
(Tab.  IX). 

The  mean  point  position  discrepancies  in  the  WGS-72 
Cartesian  coordinate  system  range  between  16.650  meters  at 
station  RM2  on  day  61  and  8.197  meters  at  station  3ZERO  on 
day  69.  The  average  of  all  the  mean  point  position  discrep¬ 
ancy  determinations  is  12.685  meters.  The  standard  errors 
of  the  mean  point  position  discrepancies  range  between  .501 
meters  at  station  BZERO  on  day  69  and  2.302  meters  at 
station  BNOB  on  day  68,  with  an  average  value  of  1.383 
meters  for  all  determinations  (Tab.  IX)  . 

E.  GEODETIC  COHPABISOHS 

The  observed  WGS-72  coordinates  of  each  station  are 
computed  by  subtracting  the  observed  discrepancies  from  the 
given  coordinates.  The  observed  W3S-72  Cartesian  coordi¬ 
nates  are  then  transformed  to  NAD-27  geodetic  coordinates 
using  a  modified  version  of  NGS’s  D035-41  transformation 
program. 

1.  Discrepancies  of  Point  Position  Determinations 

The  point  position  discrepancies  of  all  pairs  of 
determinations  are  computed  for  the  NAD-27  coordinates  (Tab. 
X) .  Azimuths  are  measured  from  north.  The  average  point 
position  discrepancy  is  -.38336  seconds  of  latitude,  -.26083 
seconds  of  longitude,  and  >9.391  meters  in  ellipsoid  height. 
The  horizontal  position  discrepancies  of  the  observed  point 
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positions  for  ail  determinations  range  between  6.498  meters 
at  station  PILOT  on  day  69  and  8.530  meters  at  station  3ZER0 
on  day  68.  The  average  distance  discrepancy  along  the 
ellipsoid  is  7.596  meters.  The  azimuths  from  the  given 
positions  to  the  observed  positions  range  between  129°  41' 
53". 04469  at  station  BZEPO  on  day  58  and  34°  43'  17". 17131 
at  station  PILOT  on  day  6S.  The  observed  positions  have  an 
average  azimuth  of  109°  25*  56". 5699  from  the  given 
positions. 

2.  Discrepancies  of  Point  Positions  for  Entire 

Recording  Sessions 

There  are  two  determinations  of  point  position 
discrepancies  on  days  61  and  68  as  a  result  of  comparing 
pairs  of  stations  in  the  computations.  An  estimate  of  the 
repeatability  of  an  entire  recording  session  can  be  calcu¬ 
lated  by  first  averaging  the  two  point  position  discrepan¬ 
cies  for  each  station  then  computing  tne  standard  error  of 
the  eight  determinations. 

The  average  discrepancies  of  entire  recording 
sessions  are  -.07980  seconds  of  latitude  and  -.25903  seconds 
of  longitude  with  standard  errors  of  .03976  and  .01764 
seconds,  respectively.  The  average  height  discrepancy  is 
8.658  meters  above  the  station  with  a  standard  error  of 
4.748  meters.  The  average  horizontal  distance  discrepancy 
is  7.538  meters  with  a  standard  error  of  .714  meters. 

3.  Relative  Position  Comparisons 

The  distances  and  azimuths  between  pairs  of  stations 
are  computed  for  both  the  given  coordinates  and  the  observed 
coordinates  by  the  inverse  computation.  The  differences  of 
two  are  computed  by  subtracting  the  given  distances  and 
azimuths  from  the  observed  values  (Tab.  XI).  Six  of  the 
seven  observed  side-lengths  are  shorter  tnan  the  distances 


TABLE  II 

Discrepancies  of  Observed  Distances  and  Azimuths 


J.D. 

61 


58 


59 


Side 

Dist .  (m . ) 

Azimuth 

JRM-RM1 

1.397 

-44* 

44". 2009 

JMR-RM2 

-0.006 

1 

o 

40" . 3464 

RM 1-RM2 

-0.231 

-10  mi 

30". 5102 

BZ  ERO-PI LOT 

-3.128 

-1". 24259 

BZERO-BNOB 

-0. 428 

-26" .09039 

BNOB-PILOT 

-0.774 

9". 49565 

BZERO-PILOT 

-4. 556 

-1". 05927 

computed  from  the  given  coordinates.  In  the  short  base  line 
tests,  the  differences  in  azimuth  are  of  the  order  of 
degrees  compared  to  differences  of  the  order  of  seconds  of 
arc  in  the  medium  base  line  tests.  In  general,  the  azimuth 
differences  decrease  as  the  side-length  differences  increase 
in  both  the  short  base  line  and  medium  base  line  tests. 


4 .  Relative  Position  Discrepancies 

Relative  positioning  is  used  to  determine  an  unknown 
position  relative  to  the  position  of  a  known  station.  The 
direct  computation  calculates  the  coordinates  of  a  position 
given  the  coordinates  of  a  known  position  and  the  distance 
and  azimuth  to  the  unknown  position.  The  distance  and 
azimuth  between  pairs  of  observed  positions  are  used  in  the 
calculation.  Stations  JMF.  and  RM 1  are  taken  as  fixed  in  the 
short  base  line  tests  and  the  positions  of  RM1  and  R M2  are 
calculated.  Stations  BZERO  and  BNDB  are  used  to  locate 
stations  ENOB  and  PILOT  in  the  medium  base  line  tests.  The 


discrepancies  of  the  calculated  positions  are  computed  in 
terms  of  latitude,  longitude,  elevation,  distance,  azimuth, 
and  relative  accuracy  between  stations  (Tab.  XII).  The 
elevation  discrepancies  are  determined  by  computing  the 
difference  of  the  height  discrepancies  between  pai~s  of 
point  position  solutions. 

Relative  position  determinations  differ  from  the 
given  positions  b'y  an  average  of  1.  452  meters  ir.  the  short 
base  line  tests  and  2.712  meters  in  in  the  medium  base  line 
tests.  The  relative  accuracies  between  stations  range 
between  1:13  for  the  side  JMR-RM1  on  day  61  and  1:17,418  for 
the  side  BNOB-PILOT  on  day  68.  The  average  relative  accu¬ 
racy  of  the  medium  base  line  tests  is  1:9,712.  Station 
PILOT  was  determined  to  relative  accuracies  of  1:8,139  and 
1:5,593  with  respect  to  station  SNOB  on  subsequent  days. 

The  data  acquired  in  the  short  base  line  tests 
provide  useful  information  in  the  absolute  position  anal¬ 
yses;  however,  this  scenario  is  not  a  practical  application 
of  relative  position  determinations.  It  is  not  justifiable 
to  measure  distances  of  the  order  of  thirty  meters  with  such 
expensive  equipment  when  the  distances  can  be  measured  accu¬ 
rately  with  steel  tape  or  other  inexpensive  equipment. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


Based  or*  the  results  of  the  medium  base  line  tests,  the 
real-time  solutions  of  the  TI4 100/3 EOSIAR  are  capable  of 
establishing  relative  positions  to  Third-order,  Class  II 
standards.  Two  receivers,  one  located  at  a  known  position 
and  the  other  separated  by  14  to  25  kilometers,  can  locate 
horizontal  control  in  real-time.  This  assumes  five  hours  of 
visibility  of  four  satellites  ana  a  means  for  the  output  of 
a  solution  based  on  the  averaging  of  real-time  solutions 
recorded  at  12-second  intervals. 

It  is  found  that  relative  positioning  using  GPS  is  more 
accurate  than  absolute  point  positioning.  The  average 
distance  discrepancy  of  th._  absolute  position  determinations 
is  7.538  meters  compared  to  relative  position  averages  of 
1.452  meters  in  the  short  base  line  tests  and  2.712  meters 
in  the  medium  base  line  tests.  The  absolute  point  position 
determinations  have  a  one  sigma  repeatability  of  .714  meters 
in  distance  discrepancy.  The  relatively  high  repeatability 
of  absolute  point  position  determinations  resulted  in  rela¬ 
tive  position  accuracies  of  1:9,712,  oq  averaye,  for 
distances  of  14  to  26  kilometers. 

Third-order,  Class  II  accuracies  can  be  obtained  with 
the  GEOSTAR  Feceiver  in  real-time.  A  typical  operational 
unit  consists  of  two,  2-man  teams;  each  team  supplied  with 
the  GEOSTAR  field  set.  One  team  sets  an  antenna  over  a 
known  position  and  the  other  team  locates  at  the  point  to  be 
surveyed.  Both  receivers  should  gather  data  simultaneously 
during  the  periods  of  satellite  visibility  (during  Phase  II 
testing).  The  teams  should  agree  to  beyin  the  averaging 
algorithm  at  a  predetermined  time  tag,  preferably  a  short 
time  after  the  uploads  to  the  satellites.  Satellite  upload 
times  can  be  requested  beforehand  from  Vandenburg  AFB. 


Additional  averaging  methods  should  be  explored.  In 
regard  to  the  real-time  solutions,  averaging  analyses  could 
be  performed  to  determine  the  amount  of  data  acquisition 
time  required  to  achieve  certain  accuracies.  An  analysis  of 
this  nature  would  be  important  for  special  military  applica¬ 
tions  where  an  accurate  position  is  required  in  less  than 
five  hours  of  receiver  operation. 

In  the  Fall  1984  the  GEOSTAB  will  undergo  hardware  and 
software  modification  before  full-scale  production.  It  is 
planned  that  the  software  will  be  modified  to  include  real¬ 
time  solutions  of  short  time-interval  observations;  these 
solutions  will  be  different  from  existing  real-time  solu¬ 
tions  based  on  a  cumulative  least  squares  method.  [fief.  22] 
Felative  position  determinations  using  the  real-time  solu¬ 
tions  should  improve  with  the  new  software  because  the  solu¬ 
tions  of  two  receivers  can  be  compared  simultaneously  with 
no  influence  from  preceding  solutions. 


APPENDIX  A 

GPS  USER-TAPE  READ-WRITE  PROGRAM 


i 


c 

C  THIS  PROGRAM  IS  A  SLIGHTLY  MODIFIED  VERSION  OF  THE 

C  EXAMPLE  INPUT  PROGRAM  CONTAINED  IN  "A  STANDARDIZED 
C  EXCHANGE  FORMAT  FOR  NAVSTAR  GPS  GEODETIC  DATA"  BY 
C  BY  V-  DAN  SCOTT  AND  J.  GARRY  PETERS,  APPLIED  RESEARCH 
C  LABORATORIES,  THE  UNIVERSITY  DF  TEXAS  AT  AUSTIN,  MARCH 
C  10,  1983.  THIS  PROGRAM  CAN  BE  EASILY  CHANGED  TO  READ 

C  AND  WRITE  ANY  DATA  TYPE  ON  GPS  USER-TAPES  SUPPLIED  BY 

C  NSWC .  TO  WRITE  ALTERNATIVE  DATA  TYPES,  CONSULT  THE 
C  ABOVE  PUBLICATION  TO  DETERMINE  THE  DATA  TYPES  THAT 
C  CORRESPOND  TO  THE  DATA  OF  INTEREST.  EXAMPLE:  TO  PRINT 

C  OUT  EPHEMERIS  INFORMATION  (DATA  TYPE  4110),  CHANGE 

C  THE  GO  TO  STATEMENT  FOLLOWING  THE  IF  STATEMENT  THAT 

C  TESTS  FOR  THIS  DATA  TYPE  BELOW  TO  TRANSFER  CONTROL 

C  TO  WRITE  STATEMENTS  .WITH  THE  PROPER  FORMAT.  THE 

C  FORMATS  OF  ALL  DATA  TYPES  ARE  GIVEN  IN  THE  ABOVE 

C  PUBLICATION.  THE  JCL  FOR  THIS  PROGRAM  IS  FOR  AN 
C  IBM  370  MODEL  3033.  THE  TAPES  ARE  RECORDED  WITH  1600 
C  BPI,  A  LOGICAL  RECORD  LENGTH  OF  30,  A  BLOCK  SIZE  OF 
C  6400,  AND  ASCII  CODED. 

C 

C 

//RAKOW  JOB  (1643,01  30)  , ' RAKOWSKY ' , CLASS=F 
//♦MAIN  CRG=NPGVM1. 1643P 
//  EXEC  FORTXCLG 
//FORT.  SYS  IN  DD  * 


c 


c 


C  C. 


[•a*; 


r 


REAL*6  TTAG,X,Y,Z,SIGT 
c 

INTEGER  ICN,  IE  EP, ICONI  F.,  NR EF, INDEX, DTYPE,  JUNK, ASTER, 

*  BLANK, SENT, PL US,CH ARE , CHAFF, CHARI, NRID 
REAL  SIGX, SIGY  ,SIGZ, SKIP 

c 

DATA  ASTER/1  *  '  /,  PLUS/*  +  ’  /,  CHA  R  E/'  2'  /,  CHARE/’  F’  /, 

*  CHAR/’  I’/,  BLANK/’  '/,SENT/'$'/ 
c 

DIMENSION  ICON  TE (IS) 

C 

100  CONTINUE 

READ  (1,10)  ICN,IREP,  ICONTP. 

10  FOFMAT(A  1,13,  1514) 

C 

120  CONTINUE 

IF  (ICN.EQ.  ASTER)  GO  TO  150 
IF  (ICN.SQ. PLUS)  GO  TO  100 
IF  (ICN.  EQ.  CHARE)  GO  TO  100 
IF  (ICN.EQ. CHAFF)  GO  TO  5  00 
IF  (ICN.  EQ.  CHARI)  GO  TO  500 
IF  (ICN. EQ. SENT)  GO  TO  500 
C 

150  CONTINUE 
C 

IF  (IEEP.  EQ.O)  IREP=  1 
NREP=0 
ICN=3LAN  K 
C 

200  CONTINUE 
C 

IND£X=0 
NREP=NREP+ 1 

IF  (IREP.  LT.O)  GO  TO  300 
IF  (NREP.  GT.  IREP)  GO  TO  100 
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GO  TO  300 

CONTINUE 

IF  (ICN.NE.  BLANK) 

GO  T 

'0 

1  20 

INDEX=INDEX* 1 

IF  (INDEX.  GT.  19)  GO  TO 

i  200 

DTYPE=ICONTR (INDEX) 

IF  (DTYPE.EQ.O)  GO 

TO 

200 

IF  (DTYPE. EQ.9 10) 

GO  TO  204 

IF  (DTYPE.  EQ.  920) 

GO  TO  ; 

201 

IF  (DTYPE.  EQ-  10  10) 

GO 

TO 

204 

IF  (DTYPE.  EQ.  10  20) 

GO 

TO 

20  5 

IF  (DTYPE.  EQ.  1030) 

GO 

TO 

204 

IF  (DTYPE.  EQ.  1040) 

GO 

TO 

20  3 

IF  (DTYPE.  EQ.  10  50) 

GO 

TO 

204 

IF  (DTYPE.  EQ.  1  1  10) 

GO 

TO 

204 

IF  (DTYPE.  EQ.  1 1  50) 

GO 

TO 

202 

IF  (DTYPE.  EQ.  20  10) 

GO 

TO 

20  1 

IF  (DTYPE. EQ. 2020) 

GO 

TO 

201 

IF  (DTYPE.  EQ.  2030) 

GO 

TO 

201 

IF  (DTYPE. EQ. 21  10) 

GO 

TO 

203 

IF  (DTYPE.  EQ.  2120) 

GO 

TO 

20  1 

IF  (DTYPE.  EQ. 21  30) 

GO 

TO 

201 

IF  (DTYPE.  EQ.  21  40) 

GO 

TO 

20  1 

IF  (DTYPE.  EQ.  22  10) 

GO 

TO 

202 

IF  (DTYPE.  EQ.  2250) 

GO 

TO 

20  1 

IF  (DTYPE.  EQ.  2270) 

GO 

TO 

201 

IF  (DTYPE.  EQ.  30  10) 

GO 

TO 

202 

IF  (DTYPE.  EQ.  32  10) 

GO 

TO 

20  1 

IF  (DTYPE.  EQ.  322 0) 

GO 

TO 

201 

IF  (DTYPE.  EQ.  3250) 

GO 

TO 

20  1 

IF  (DTYPE.  EQ.  40  10) 

GO 

TO 

201 

IF  (DTYPE.  EQ. 4050) 

GO 

TO 

20  1 

IF  {DTYPE.  EQ.  41  10)  GC  TO  204 
IF  (DTYPE.EQ.  4120)  GO  TC  20  2 
IF  (DTYPE.  EQ.  42  10)  GO  TO  203 
IF  (DTYPE.  EQ.  4250)  GO  TO  202 
IF  (DTYPE.  EQ.  4270)  GO  TO  202 
IF  (DTYPE.EQ.  50  10)  GO  TO  20  1 
IF  (DTYPE.EQ. 60  10)  GO  TO  6010 
IF  (DTYPE.  EQ. 6020)  GO  TO  202 
IF  (DTYPE.  EQ.61  00)  GO  TO  201 
IF  (DTYPE.EQ. 61  10)  GO  TO  20  3 
IF  (DTYPE.EQ. 70  10)  GO  TO  201 
IF  (DTYPE.  EQ.  71  10)  GO  TO  202 
IF  (DTYPE.  EQ. 71  20)  GO  TO  201 
IF  (DTYPE.  EQ.  72  10)  GO  TO  20  1 
IF  (DTYPE.EQ. 7220)  GO  TO  201 
IF  (DTYPE.  EQ.  7230)  GO  TO  201 
IF  (DTYPE.  EQ.  72  40)  GO  TO  202 
IF  (DTYPE.  EQ.  7250)  GO  TO  202 
IF  (DTYPE.  EQ.  7260)  GO  TO  203 
WRITE  (4,31) 

FORMAT  ('  DATATYPE  ERROR') 


201  CONTINUE 


READ  (1,11)  ICN 
FORMAT  (A  1) 

IF  ((ICN.  ME. BLANK)  .AND.  (ICN.  NS.  PLUS)  )  GO  TO  120 
GO  TO  300 


202  CONTINUE 


READ  (1,11)  ICN 

IF  ( (ICN.  NE. BLANK)  .AND.  (ICN. NE.  PLUS)  )  GOTO  120 
READ  (1,11)  NCN 
GO  TO  300 


203  CONTINUE 
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BEAD  (1,11)  ICN 

IF  ((ICN. NE. BLANK)  .AND.  (ICN. NE. PLUS)  )  GO  TO  120 
EE  AD  (1,11)  NCN 
READ (1,11)  NCN 
GO  TO  300 

204  CONTINUE 

BEAD  (1,11)  ICN 

IF  ( (ICN.  NE.  BLANK)  .AND.  (ICN.  NE.  PL'JS)  )  GO  TO  120 

READ (1,11)  NCN 

READ  (1,11)  NCN 

READ (1,11)  NCN 

GO  TO  300 

205  CONTINUE 

READ  (1,11)  ICN 

IF  ((ICN.  NE. BLANK)  .AND.  (ICN. NE. PLUS)  )  GO  TO  120 

READ  (1,11)  NCN 

READ  (1,1  1)  NCN 

R  EAD (1,11)  NCN 

READ  (1,11)  NCN 

GO  TO  300 

6010  CONTINUE 

READ { 1 , 1 2)  ICN,NRID,TTAG,X,Y 

IF  ( (ICN.  NE.  BLANK)  .  AND.  (ICN.  NE.  PLUS)  )  GO  TO  120 
READ (1,13)  NCN , Z, SIGT, SIGX , S I G Y , S IGZ , S KIP 
WRITE  (4, 14)  TTAG,X, Y,Z 
WRITE (4 , 15)  SIGX,SIGY,SIGZ 

READ  FORMAIS 

12  FORMAT (A  1 ,12,8 X,D22. 15,D22.15,D17.10) 

13  FORMAT  (A1,2D17.  10,4E10.3) 


C  WRITS  FOF MATS 
C 

14  FORMAT (D22. 15, D22. 15,D17.10,D17.10) 

15  FORMAT (31X,D17.10,4X,E10.3,4X,E10.3) 

C 

GO  TO  300 
500  CONTINUE 
STOP 
END 
c 
c 

//GO. FT01F001  DD  UNIT=3400-3 , VOL=SER=my tape , 

//  DISP=  (OLD, PASS) , LABEL= (2,NL,,IN) , 

//  DCB  = (RECFM=FB,LRECL=80,BLKSIZE=6400,DEN=3,0PTCD 
//GO. FT04F001  DD  DSN =MSS. SI  643 . FILENAME, UNIT= 333 0  V , 
//  DISP= (NEW, CATLG, DELETE)  ,  MSVGP=P0B4B, 

//  BCB= (RECFM=FB, LR ECL=80 , BLKSIZE=6 4 00 ) 


APPENDIX  3 

PROGRAM  TO  COMPOTE  REAL-TIME  SOLOTION  DISCREPANCIES 


£*******************  ********  ***************************  **  * 

c 

C  THIS  PROGRAM  COMPUTES  THE  DIFFERENCES  OF  THE  X,  Y, 

C  AND  Z  COORDINATES,  THE  TOTAL  POINT  POSITION  DISCREP- 
C  ANCIES,  AND  THE  SIDE-LENGTH  DI5 CREPANCIES  FOR  TWO 
C  STATIONS  AT  COMMON  GPS  TIMES.  THIS  PROGRAM,  IN  EFFECT, 
C  WRITES  THE  ABOVE  DATA  BEGINNING  WITH  THE  FIRST  REAL- 
C  TIME  SOLUTIONS  OF  EACH  TAPE,  THEN  SPOOLS  THE  TAPE  OF 
C  THE  EARLIER  RECORDING  SESSION  TO  CORRESPOND  TO  THE 

C  GPS  TIME  OF  THE  TAPE  THAT  STARTED  RECORDING  AT  A  LATER 

C  TIME.  A  CONVERSION  OF  THE  GIVEN  GEODETIC  COORDINATES 
C  TO  WGS-72  CARTESIAN  COORDINATES  MUST  FIRST  BE  DONE. 

C  THESE  VALUES  ALONG  WITH  COMPUTED  SIDE-LENGTH  DISTANCE 
C  ARE  ENTERED  BELOW  FOLLOWING  THE  TYPE  DEFINITIONS.  THE 
C  OUTPUT  IS  IN  THE  FORM:  ’GIVEN  VALUE  -  REAL-TIME 
C  SOLUTION  VALUE’. 

C 

C* *************  ********  *********************************** 

c 

//RAKOW  JOB  (1643,0130)  ,  '  RAKOWSKY*  ,CLASS=B 
//♦MAIN  ORG=NPGVM1. 1643P 
//  EXEC  FORTXCLG 
//FORT.  SYSIN  DD  * 
c 

RE AL*8  TTAG1 ,TTAG2, XI ,X2, Y1, Y2,Z1 ,  Z2, DELX, DELY,DELZ , 

*  DIST, D  VA  LUE , SIGX1 ,SIGX2, SI3 Y 1 , SIGY  2, SI3Z 1 , SIGZ2 , 

*  DELDST,X 1GIVN,X2GIVN, Y1GIVN, Y23I VN , Z 1GI V N , Z2GIVN, 

*  XI  DIF, X2DIF, Y1DIF, Y2DIF,Z 1DIF,Z2DIF, OFSET1 ,OFSET2 
C 

D VALUE  =  25  493.  188 


» 


» 


I 


I 
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X1GIVN=-740577.548 
rlGI7N=-5456773.662 
Z1GIVN =3207660 .081 
X2GI7F=-740336 .50 1 
Y2GIVN=-54697 17.020 
Z2GI VN=3 185698. 43 2 

READ (1 ,1 0, END=500)  IT AG  1 , X 1 , Y  1 , Z 1 
READ  (1*11)  SIGX1,SIGY1  , SIG  Z 1 
READ  (2,10, END= 500)  TTAG2, X2, Y2, Z2 
RE  AD  (2, 11)  SIGX2,SIGY2,SIGZ2 

IF  (TTAG1  .LT.  TTAG2)  GO  TO  3 

2  CONTINUE 

READ (2,10)  TTAG2, X2 , Y2 , Z  2 
READ  (2,11)  SIGX2,SIGY2,SIGZ2 
IF  (TTAG2.GE.TTAG1)  GO  TO  30 
GO  TO  2 

3  CONTINUE 

READ  (1,10)  ITA G1,X1,Y1,Z1 
R  EAD (1,11)  SIGX1,SIGY1 ,SI3Z1 
IF  (TTAG1.GE.TTAG2)  GO  TO  30 
GO  TO  3 

30  CONTINUE 

IF  (TTAG1  .  GT.  TTAG2)  GO  TO  2 

IF  (TTAG2.GT.  TTAG1)  GO  TO  3 

DEIX=X1-X2 

DEIY=Y1-Y2 

DELZ=Z1-Z2 

DIST  =  DS2RT (DEL  X**2*DEL  Y*  *2  +  DELZ**2 ) 
DELDOr:=DIST-DVALUE 


X 1  DI  F=  XI  -  X  1 G I V  N 
Y1DIF=Y1-Y1GIVN 
Z1DTT'=Z1-Z1GIVN 
X2DIF=X2-X2GIVN 
Y2DIF=Y2-Y2GIVN 
Z2DIF=Z2-Z2GIVN 

0FSET1 =DSQRT  (X1DIF**2  +  Y1DIF**2+Z1DIF**2) 
OFSET2=DS2RT (X2DIF**2*Y2DIF**2*Z2DIF**2) 
C 

WF.ITE  (4,  42)  TTAGl ,  TTAG  2 
WRITE  (4,43)  XI DIF,  Y 1  DI  F, ZlDIF 
PR  ITE  (4,41)  0FSET1 
WRITE  (4,43)  X2DIF, Y2DIF, Z2DIF 
WRITE  (4,41)  0FSET2 
WRITE  (4,41)  DFIDST 
C 

READ  (1,10,  END=500)  TTAG1 , X 1 ,  Y  1 ,  Z 1 
E  EAD (1,11)  SIGX1,SIGY1,SIGZ1 
READ (2,10, END=  500)  TTAG2 ,X2,Y2,Z2 
READ (2, 11)  SIGX2,SIGY2,SIGZ2 
GO  TO  30 
C 

C  READ  FORMATS 
C 

10  FORMAT  (D22.  15,  D22.  1 5,  D  1 7 .  10,  D  1 7.  1  0) 

1 1  FORM  AT (31X,D17 . 10,4X,E10.3,4X,E10. 3) 

C 

C  WRITE  FORMATS 
C 

41  FORMAT  (D22.  15) 

42  FORMAT  (2D22.  15) 

43  FORMAT  (3D22.  15) 

C 

500  CONTINUE 
STOP 
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//30. FT01F001  DD  DSN=MSS. SI  643 . MYIAPE, 

//  DISP= (OLD, KEEP) , MSVGP= PUB45 

//30 . FT02F001  DD  DSN=MSS. SI  6 43. FILENAME, 

//  DISP= (OLD, KEEP) ,MSVGP=PD34B 

/* 30 . FTO  4F00 1  DD  DSN =MSS . SI  643 . CHE- K , UNI T= 33 50 , 

/*  DISP=  (NSW,  KEEP)  ,VOL=SER=MVS034,SPACE=  (CYL,  (4,4)) 

/*  DC  E=  (F.SCFtf  =  FB,LRECL=80,BLKSIZE  =  6400) 

/* 

//30. FT04F0C1  DD  DS N=MSS . S 1 6 43 . filename, UNIT=333 OV , 

//  DISP= (NEW, CATLG, DELETE) , MSVGP=PU34 B, 

//  DCB= (RECFM=FB,LRECL=80,BLKSIZE=64Q0) 
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cdu 

DMA 
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GEOSTAR 

Gil! 

GPS 

Mbps 

MCS 

NAD- 2 7 

NAVSTAR 

NSWC 

N SWC/DL 

PPS 

PRN 

SPS 


Age  of  Data  (Clock) 

Age  of  Data  (Ephemeris) 
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Defense  Mapping  Agency 

Defense  Mapping  Agency  Hydro. /Topo.  Center 

Geodetic  Satellite  Tracking  and  Ranging 

Greenwich  Mean  Time 

Global  Positioning  System 

Million  bits  per  second 

GPS  Master  Control  Station 

North  American  Datum  of  1927 

Navigation  Satellite  Tracking  and  Ranging 

Naval  Surface  Weapons  Center 
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Precise  Positioning  Service 

Pseudo-Random  Noise 

Standard  Positioning  Service 
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